Effect of ionization/recombination processes on the electrical 
interactions between positively charged particles in highly 

collisional plasmas 



M. Chaudhuri, S. A. Khrapak, and G. E. Morfffl 

Max-Planck Institut fur extraterrestrische Physik, D- 857 4-1 G arching, Germany 

(Dated: October 20, 2010) 

The effect of ionization and recombination processes on the electrical interactions 
between a pair of small charged particles in highly collisional plasmas is discussed. 
In particular, it is shown that these processes suppress the long-range attraction 
between positively charged particles. The condition corresponding to the vanishing 
of attraction is derived. The role of the effect for conditions of existing experiments 
is estimated. 

The shape of the interparticle interaction potential is a very important factor which deter- 
mines various physical phenomena in complex (dusty) plasmas (as well as in other interacting 
particle systems). Different mechanisms of attraction and repulsion between the charged par- 
ticles immersed in plasmas have been discussed in the literature." At short distances electri- 
cal repulsion between like-charged particles dominates. In isotropic conditions this repulsion 
can be approximated by the Debye-Hiickel (Yukawa) potential U(r) = {Q 2 /r) exp(— r/A), 
where Q is the particle charge and A is the plasma screening length. At longer distances 
various effects can contribute to the interparticle interaction. Among these is the so-called 
"ion shadowing" interaction,-^ which can generate an attractive branch of the interaction 
potential.- 1 ^ A similar effect known as "neutral shadowing" can lead to additional attractive 
or repulsive contribution to the interparticle interaction potential, depending on the sign of 
the temperature difference between the particle surface and surrounding gas.— ^ Even the 
potential of electrical interaction between the charged particles can have complicated struc- 
tures at large interparticle distances. In general it is not screened exponentially, but has an 
inverse power-law long-range asymptote as long as plasma production and loss processes in 
the particles vicinity are neglected.- 1 ^ The potential decays as cx 1/r 2 in collisionless plas- 
masi3aj an( j as oc 1/r in highly collisional plasmas.— iliLIS Moreover, electrical interaction 
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can exhibit long-range repulsion or attraction, depending on the sign of the particle charge. 

As long as collective effects relevant to dense dust clouds are neglected, two negatively 
charged particles repel each other electrically at any distance. However, two positively 
charged particles can attract each other at long distances. The possibility of attraction 
between individual positively charged emitting particles was pointed out by Delzanno et al. 
who performed a simulation study of the particle charging in a collisionless plasma, taking 
into account the effect of thermionic emission from the particle surface.— They constructed 
an analytical theory to describe this effect in collisionless plasmas.— Later on, an analytical 
approach to calculate the attractive part of the potential between positively charged emitting 
particles was developed for the highly collisional plasma regime. 19 The physical mechanism of 



this attraction was discussed in Ref. 
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and it was suggested there that attraction can explain 
the formation of ordered structures in the "dusty combustion experiment". 20 However, an 
important effect - ionization and recombination in the background plasma - was neglected 
in this theoretical study. 

As has been recently demonstrated, the potential of electrical interaction between the 
particles is sensitive to the strength and exact nature of the plasma production and loss 
processes.—* 2 ^ So far, only the case of negatively charged particles has been considered in 
literature. The purpose of this Brief Communication is to estimate the effect of ionization 
and recombination processes on the attractive branch of the interaction potential between 
a pair of small positively charged particles in plasmas. In particular, we find that ionization 
suppresses the attractive part of the potential and, when the ionization rate becomes suffi- 
ciently high, the attractive branch disappears. We identify the condition when this happens. 
We also estimate the importance of ionization/recombination processes for the conditions of 
the "dusty combustion experiment" , and demonstrate that here they play only a negligible 
role. 

We assume, below, that electron impact ionization is responsible for the plasma pro- 
duction, while plasma losses are due to the electron-ion volume recombination.—^ 2 - In the 
limit of highly collisional plasma (li< e ) <C Ad), both the ion and electron components are 
mobility controlled and can be described by hydrodynamic equations. Here £^ is the 
mean free path of ions (electrons), Ad = (A^ 2 + A^ 2 ) -1 ^ 2 is the linearized Debye radius, 



Aoi(e) = ^/Ti( e )/4:7in e 2 is the ion (electron) Debye radius, is the ion (electron) tem- 
perature, and rii ~ n e ~ n is the unperturbed plasma density. Due to the presence of 
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fast thermalization processes in highly collisional plasma, we can assume that the electron 
temperature (of both emitted and background electrons) is uniform. Finally, we consider the 
potential distribution around an individual particle and the resulting interaction between a 
pair of particles, thus neglecting any collective effects including plasma losses on the particle 
component. 

The distribution of electrical potential around an individual particle can be calculated 
combining the continuity and momentum equations for electrons and ions in the hydrody- 
namic approximation with the Poisson equation. The result is^ 1 ^ 



where 



and 



{Q+/r) exp(-rk + ) + (Q-/r) exp(— rfc_), (1) 



Q [k% - kl - (eJo/QDi) (1 - A/Avj 
Q± = T fc 2 _ k 2_ — 1 



2k 2 ± = kl + v I /D i ± sj {kl + vt/DiY - Auj {klJD t + k 2 Di /D e ). (3) 

Here A(e) = v T i( e )ti(e), v T i( e ) = ^T iie) /m^ e) , m i(e) , and fc Dl -( e ) = A^* (e) are the diffusion 
coefficient, thermal velocity, mass, and inverse Debye radius for ions (electrons) respectively. 



The inverse linearized Debye radius is k^> = \/^De + ^m- the stationary state the electron 
and ion fluxes collected by the particle are equal to each other, J, = J e = Jo- The total 
electron flux J e includes both the contribution from the electrons absorbed on the particle 
surface and those emitted from the surface. We assume here that electron emission is quite 
significant, so that the particle charge is positive, Q > 0. The ionization rate vi is related 
to the recombination constant via vj = /3no- 

Equation ([T|) demonstrates that in the considered case the potential is screened exponen- 
tially, but unlike in the Debye-Huckel theory the potential contains the superposition of the 
two exponentials with different inverse screening lengths (screening parameters) k + and k-. 
Both of these screening parameters depend on the strength of plasma production through 
the ionization frequency Uj. The effective charges Q + and Q_ both depend on the strength 
of plasma production and plasma fluxes collected by the particle. The short range part of 
the potential is determined by the first term in Eq. ([1]) with the effective charge Q + and 
the larger screening parameter k + . It can be easily shown that Q + > and therefore this 
part is always repulsive. On the other hand, the long range asymptote of the potential is 
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determined by the smaller screening parameter with the effective charge Q_. The long- 
range attraction operates when Q_ < 0. From Eq. (j2J) we immediately get an approximate 
condition for the existence of attraction. It can be written as 

Jo>MQ/e). (4) 

It is obvious from the obtained inequality that attraction between two positively charged 
particles exists when the ionization strength is low enough. Otherwise, the potential of 
electrical interaction is purely repulsive. 

Note that if we neglect completely the ionization and recombination processes (uj = 
0,(3 = 0) the potential (JTj) reduces to that obtained previously [Eq. (8) of Ref. |l9|. The 
effect of ionization is small when > uj/Di. In this case the smaller screening parameter 
which determines the global screening is k 2 _ ~ k^iyi / k^D^l + (k^Di/k^De)}. Usually 
(k^Di/k^De) <C 1 and, therefore, fc_ ~ ^DeV 'vi/k^Di. This yields a lengths scale R ~ 
^De^/k^Di/vj > Aoe above which ionization and recombination in plasma should be taken 
into account. For r < R the model neglecting ionization and recombination processes is 
applicable. 



As suggested in Ref. [19|, electrical attraction between positively charged emitting particles 
can explain the formation of ordered structures in the "dusty combustion experiment" by 
Fortov et al.—^^ In this experiment the dust particles (Ce02 grains of radius a ~ 0.4 /im) 
were injected into a laminar air spray at atmospheric pressure and temperature T ~ 1700 — 
2200 K created by a two-flame Meeker burner. The complex (dusty) plasma constituents 
were air, electrons, Na + ions and Ce02 particles, all in thermal equilibrium. The grains were 
charged positively by emitting thermal electrons up to Q ~ 10 2 e. At T ~ 1700 K the particle 
component formed a short range ordered structure, with the pronounced first maximum in 
the pair correlation function. This phenomena can be attributed to the presence of the weak 
attractive part of the interparticle interaction potential, with the minimum of the potential 
corresponding approximately to the position of the first maximum in the pair correlation 
function.^ However, the effects of ionization and recombination were not taken into account 
in Ref. U 



Let us therefore estimate their importance for the considered experimental 
conditions. 

Let us first estimate the important plasma parameters. Some of them were already 
estimated in Ref. |l9|. According to these estimates for T ~ 1700 K and n e ~ 7 x 10 10 
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cm~ 3 we have £i ~ 0.06 /im, Ad ~ 10 fim, and Vt% ~ 8 x 10 4 cm/s. This yields Di ~ 
0.5 cm 2 /s for the ion diffusion coefficient. The next step is to estimate the ionization 
rate vj for the condition of the experiment. We use the Saha equation,— (n e rii/N) = 
2(gi/g a )(2iim e T/h 2 ) 3 / 2 exp(— I/T), to determine the density of Na atoms. Using gi = 1, 
g a = 2, and I ~ 5.1 eV (here / denotes the ionization potential of Na atoms) we get 
iV ~ 2 x 10 16 cm -3 . Then, an approximate expression for the electron impact ionization 



frequency,-^ v I ~ + 2T) exp(-//T), yields vi ~ 4 x 10~ 7 s _1 , where d ~ 

10~ 16 cm 2 /eV is used.— Finally, for the ion flux collected by the particle we get 1 ^ Jo — 
47ran -DjZexp(— z) ~ 4 x 10 6 s _1 , where z = Qe/aT ~ 2.5 is the dimensionless particle 
surface potential. Looking at these numbers we immediately see that fc D 3> fi/D^ and 
Jo ^ vi{Q/e). This implies that for this particular experiment, the effect of ionization 
and recombination plays a negligible role and the conclusions of Ref. |l9| remain unaffected. 
We would like to point out that, in the experiment discussed above, plasma losses on the 
particle component can play a considerable role due to the large plasma fluxes directed to 
the particles and non-vanishing particle concentration. This effect along with potentially 
important collective effects is, however, beyond the scope of this Brief Communication and 
will be considered separately. 

To summarize, ionization and recombination processes constitute an important factor 
which affects the distribution of electrical potential around an individual particle in plasmas 
and therefore is important for electrical interaction between the particles. We have shown 
here that ionization/recombination processes can suppress the attractive branch of the in- 
teraction potential between a pair of positively charged particles, and derived the condition 
when this attraction vanishes. According to our estimates, ionization/recombination effects 
are negligible for the conditions of the "dusty combustion experiment" . However, we cannot 
exclude that in other regimes they can play a dominant role. 
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